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Magnetostatic (dipolar) interactions between nanoparticles promise to open new ways to design 
nanocrystalline magnetic materials and devices if the collective magnetic properties can be controlled at the 
nanoparticle level. Magnetic dipolar interactions are sufficiently strong to sustain magnetic order at 
ambient temperature in assemblies of closely- spaced nanoparticles with magnetic moments of > 100 \i B . 
Here we use electron holography with sub-particle resolution to reveal the correlation between particle 
arrangement and magnetic order in self- assembled ID and quasi-2D arrangements of 15 nm cobalt 
nanoparticles. In the initial states, we observe dipolar ferromagnetism, antiferromagnetism and local flux 
closure, depending on the particle arrangement. Surprisingly, after magnetic saturation, measurements and 
numerical simulations show that overall ferromagnetic order exists in the present nanoparticle assemblies 
even when their arrangement is completely disordered. Such direct quantification of the correlation between 
topological and magnetic order is essential for the technological exploitation of magnetic quasi-2D 
nanoparticle assemblies. 
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The assembly of magnetic nanoparticles into larger structures is of intense interest as an approach to engineer 
new materials and devices, such as strong magnets for high-end electronics 1 , logic gates 2 , thermally stable 
magnetic memories 3 and soft magnetic materials for magnetic sensors 4 ' 5 . It also provides new opportunities 
for shaping and diversifying components in devices, by making use of magnetic interactions between particles to 
build elaborate anisotropic structures such as chains and dendrites 6 " 8 . However, in order to reveal the full 
potential of these materials, quantification of the interplay between nanoparticle arrangements and magnetic 
order is essential. 

Dipolar ferromagnetism ("magnetostatic alignment of all magnetic moments") was first predicted by Luttinger 
and Tisza for a face- centered cubic lattice of point dipoles 9 . It was subsequently explored both theoretically and 
experimentally for nanoparticle assemblies, including particle chains 10 , ordered 2D arrays 11 " 15 , and 3D lattices 16 . 
Dipolar glasses in random 3D arrangements 13,17 " 19 and dipolar antiferromagnetism in 2D arrays were also 
investigated 20 . However, local quantitative observations of dipolar magnetism in nanomaterials are lacking 
due to stringent experimental requirements. Consequently, the physical picture behind dipolar magnetism in 
nanoparticulate materials has developed with limited experimental support. 

In conventional materials, dipolar magnetism is generally considered negligible, as the classical magnetostatic 
interactions between neighboring atoms with magnetic moments of only a few |i B > a few A apart, are too weak to 
explain magnetic ordering. Without electronic exchange coupling between atoms, solid-state magnetism ceases at 
non-vanishing temperatures. However, in "meta-materials" where atoms are replaced by nanoparticles, the 
"atomic" description of magnetism (i.e., strong exchange coupling and negligible dipolar interaction) breaks 
down as dipolar interactions may overwhelm inter-particle exchange coupling leading to dipolar magnetic 
ordering at 100-1000 K 410 " 23 . 

An assembly of nanoparticles can be considered as a meta-material whose structure can range from crystalline 
when the underlying nanoparticle lattice is regular to amorphous when the particles are distributed randomly. 
The energy scale of dipolar interactions between magnetic nanoparticles varies with |i 0 |i 2 /47rr 3 and is —2500 K for 
two 15 nm diameter Co particles with moments |i of ~2* 10 5 |i B and a center-to-center separation of r = 20 nm. 
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This value is comparable to the exchange coupling energy of 1388 K 
for Co, and about four orders of magnitude larger than the dipolar 
coupling between two atomic spins at typical interatomic distances in 
a crystal. 

Just as magnetic order and properties depend on the arrangement 
of atoms in a conventional material, so are dipolar magnetic order 
and properties expected to emerge from the geometry of an assembly 
of interacting nanoparticles in a nanoparticulate material; however, 
in such materials local structure-property relationships have not 
previously been quantified experimentally in detail. We have inves- 
tigated dipolar magnetism at high spatial resolution by preparing a 
set of low-dimensional structures of self-assembled 15 nm cobalt 
particles, and using transmission electron microscopy and electron 
holography as primary imaging techniques. As a result of the aniso- 
tropic nature of the particles' magnetic dipole fields, elongated struc- 
tures such as chains may form 6 " 8,24 and local magnetic ordering can 
depend sensitively on the geometric arrangements of neighbouring 
particles: For instance, ferro- and antiferromagnetic dipolar ordering 
are predicted to occur for triangular and square particle lattices, 
respectively 20 (see Fig. 1). Electron holography 25,26 provides quant- 
itative access to nanoscale magnetic field distributions with a spatial 
resolution well below the particle size. It is currently the only tech- 
nique that can resolve non-invasively dipolar magnetism at this high 
resolution and sensitivity. 

Here, we present the first direct observations of the local magnetic 
dipolar ordering in low-dimensional nanoparticle structures, and we 
correlate this to the exact geometric arrangement of the particles and 
to the magnetic history of the assemblies. The results clearly show the 
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Figure 1 | Simulation and schematic diagrams illustrating magnetic 
dipolar fields and interactions between spherical nanoparticles. 

(a) Magnetic field cross-section of the end of a chain of 10 uniformly 
magnetized spheres (red outline). The strongest field occurs at its ends. The 
fringing field gradient may guide additional particles to attach 
preferentially at the ends of the chain in a ferromagnetic alignment (top 
particle with black outline) . At the side of the chain the strong field between 
two particles favors ferromagnetic alignment if the added particle (black 
outline) is located within a relatively small region (labeled FM in a), while 
antiferromagnetic (AFM) alignment will occur otherwise (particle with 
white outline). For a freely moving particle the close-packed lattice with 
ferromagnetic alignment is energetically favored, (b) Schematic 
illustrations of triangular (close-packed) and square lattices leading to FM 
and AFM dipole ordering, respectively. 



existence of dipolar magnetic ordering, which interestingly takes 
different forms depending on particle arrangement and magnetic 
history. Far from being "weak", dipolar magnetism in low-dimen- 
sional nanoparticle structures reveals a tendency to survive major 
degrees of structural disorder in the arrangement of the nanoparti- 
cles, to the point that we envision the existence of amorphous meta- 
materials with ferromagnetic character. 

Results 

We prepared assemblies of 15 ± 2 nm s-Co particles with an oleic 
acid coating by evaporation of a colloid 8,27 on a carbon substrate 
(with no external field applied). Figures 2 and 3 show transmission 
electron microscopy (TEM) bright-field images of particle assem- 
blies. For chains wider than 1 particle across, particles are typically 
assembled into triangular (close-packed) lattices, although square 
lattice arrangements are also occasionally seen. The particles are 
separated by at least 2 nm (the size of the surfactant layer) providing 
full exchange decoupling. 

We applied off- axis electron holography 25,26 to map directly, and 
non-invasively, the projected magnetic field, #i n - p iane> of the elon- 
gated nanoparticle assemblies with a nominal spatial resolution of 
6.3 nm (fringe spacing is 2.1 nm). Figure 2 shows color-coded holo- 
graphic field maps acquired at remanence before and after applying 
an off-plane field of ± 2 T, revealing the magnetic moment topo- 
graphy of the chains. In a single-particle- wide chain (i), simple fer- 
romagnetic order exists already in the initial state (where chains may 




Figure 2 | Bright-field (BF) TEM images and holographic magnetic field 
maps of chains of ~ 15 nm Co particles that are (i) one, (ii) three and (iii) 
multiple nanoparticles in width. Initial remanent magnetic states (IS) and 
magnetic states after application of a field of ±2 T are shown. The white 
contours and the colors (described using the color wheel) provide a 
measure of the magnitude and direction of the projected magnetic field in 
the plane of the particles with a spatial resolution of ~6 nm. Each white 
contour represents a phase variation of 7i/64, corresponding to 1/64 of a 
flux quantum 4> 0 = 2.07* 10" 15 T*m 2 enclosed between two adjacent 
contours. 
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Figure 3 | (a) Bright-field TEM images and (b) holographic magnetic 
field maps obtained from —15 nm Co particle chains of increasing width 
numbered I— VII at remanence after magnetic saturation. CN= 

coordination number (= average number of nearest neighbors), a g = 
geometric disorder, a m = magnetic disorder (angle measured in radians), 
M = relative net magnetic moment of chain. Phase contour spacing as in 
Fig. 2. 
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have been exposed only to the weak stray field, ~6 mT, of the micro- 
scope column while inserting the sample). Chains wider than one 
particle show more complicated initial states. In a 3 -particle- wide 
chain (ii), ferromagnetic alignment is observed only for two of the 
three strands while antiferromagnetic alignment is found between 
these two strands and the third one. The antiferromagnetic align- 
ment seems to be related to cases of square lattice arrangement 
between some of the particles. In a 6 -particle- wide chain (iii), several 
flux-closed regions with domain widths of up to 5 particles are 
observed in the initial state. No out-of-plane cores, typical for regular 
vortices in continuous thin films, are found as a result of the absence 
of inter-particle exchange interactions. At remanence after magnetic 
saturation in a field of ± 2 T, overall dipolar ferromagnetic order is 
observed along the direction of chains i-iii, with only local short- 
range variations ("ripples") of the moment orientations in chains ii 
and iii. 

In order to quantify dipolar ferromagnetic order in the nanopar- 
ticle structures we examine a set of chains (I-VII, Fig. 3) with widths 
of 1-8 particles after applying a field of 2 T. The average number of 
nearest neighbour particles (coordination number, CN) of each par- 
ticle ranges from 2 for the single-particle chain (I) to 5.65 for the 
wider chains corresponding to quasi-2D particle arrangements (CN 
= 6 for infinite close-packed 2D lattices). 

From the projected field maps (Fig. 3b) we estimate the magnitude 
and orientation of the magnetic moment of each individual particle 
(Fig. 4a shows this for chain IV) and correlate it with the geometrical 
arrangement of the particles. For each pair of particles in each chain 
we measure the spatial separation r between their centers and the 
angular difference A0 between their moments. Figure 4b shows the 
distribution (r,A6) for chain IV. We found for each chain the dimen- 
sionless geometric order parameter ratio a g = a nn /d nn , where d nn is 
the average nearest neighbor distance and a nn is the spread in d nn . 
We also found the degree of local magnetic disorder rj m from the 
spread in A0 between nearest neighbors, and the overall relative net 
moment M per chain as the actual vector sum of all moments in a 
chain divided by the nominal sum of all moments. 



Figure 4 | Results obtained from measurement of chain IV in Fig. 3. 

(a) Estimated magnetic moment distribution, (b) radial- angular (r,A6) 
distribution function, where r is distance between two magnetic moments 
(particle centers) and AG is the angular difference in their moment 
orientations, measured for all pairs of moments in the chain. The thin 
green circular lines superimposed on (b) indicate the expected locations of 
the first few neighbouring peaks for a 2D close-packed lattice with a 
measured nearest neighbour distance d nn of 19.8 nm. 

As given in Fig. 3, we find that M indicates almost perfect fer- 
romagnetic order (M = 0.97) for the single particle chain; also, a m is 
small. For wider chains, M is ~ 10% lower and a m is about two times 
larger, but a m and M reveal that the magnetic order is strongly 
ferromagnetic on both short- and long-range scales. All chains have 
a similar geometric disorder rj g of around 0.17 ± 0.04. The spread in 
particle size (relative to mean particle size) gives values similar to 
those of a g for all chains. This suggests that the geometric disorder 
may be driven primarily by the size distribution, and that the packing 
of the particles on the substrate is almost as dense as possible for the 
given particle size distribution. 

In order to simulate magnetic order in the chains, we used a 
computational framework based on an adaptation of the Landau- 
Lifshitz-Gilbert equations that describe semi-classically the dyna- 
mics of interacting magnetic moments inside materials 28 . Starting 
from the measured coordinates and radii of the particles, and assum- 
ing an initial in -plane alignment that mimics the external field 
applied to saturate the structures, we calculated the expected ground 
states upon removal of the field for particles with zero magnetic 
anisotropy. Despite edge effects and a few minor discrepancies 
between simulations and observations, for example in the exact posi- 
tion of the vortex in chain IV, the simulations confirm magnetic 
ordering as being long-range dipolar ferromagnetic. 

We also simulated relaxation from an initial random distribution 
of moments, obtaining results similar to the initial magnetic states, 
with shorter range magnetic order that includes flux- closed regions 
that are a few particles in width for the chains wider than 1 particle. In 



SCIENTIFIC REPORTS | 3 : 1234 | DOI: 1 0.1 038/srep01 234 



3 



®®®® 
®®®® 

®®®® 
®®®® 

®®®® 
®®®® 



®®®® 
®®®® 

®®®® 
®®®® 

®®_®_®_ 



®®®® 
®®®® 

e®®® 

a g =0 



®e©e 
®®0 0 

00 ®® 
®®<s> ® 

®®®®^ 

®®®<£ 

0® <S>® 

<D<2>0®_ 

®(S>®® 

®0® ® 
®®®0 
®€>®<D 

(2)®®® 
®0®(£> 

®®® ® 
®®0® 

®®0® 
®® ®® 

j2)<D<a>® 
®®®® 

® ®0® 

©®e® 



a g = 0.2 



<5>®(5 <5> 

@® ® 0 
€)€> <i> ® 

® <S>(5®_ 
®®00 

®e®0 

® ®€>0 

®0®® 

®®® «L 
<f>®®<* 

®*> ® ® 
<5® ® ® 

>®®® 
® ®®® 

® ®®(S) 

0® e® 



a g = 0.4 




Figure 5 | Simulations of magnetic order in chains with increasing lattice disorder, (a) 4 X 18 particle chain with perfect triangular lattice (a g = 0), 
(b, c) 4 X 18 particle chain with disordered particle arrangement (a g = 0.2-0.4), (d) plot of the magnetic disorder a m vs. geometric disorder a g , 
obtained from a set of 35 simulations with varying a g , (e) plot of the normalized net magnetic moment M as a function of geometric disorder a g . 
Color-code as shown in Fig. 4. 



these cases, many degenerate magnetic states exist. This correspond- 
ence between simulated and observed states together with the pres- 
ence of complex initial magnetic states suggest that the assembly into 
chains is not entirely magnetically driven. 

The simulations were used also to assess structures that were not 
available experimentally by artificially increasing the degree of geo- 
metrical (lattice) disorder for a finite lattice of dipoles by adding a 
random component to the position and size of each nanoparticle. An 
example of a chain with dimensions of 4 X 18 nanoparticles (CN = 
4.85), is shown in Figs. 5a-c, where the initially regular triangular 
lattice (rjg = 0) is distorted gradually up to a g = 0.4 such that the 
particle arrangement is indistinguishable from an amorphous struc- 
ture, i.e. order at nearest neighbors is lost. Even in this extreme case, 
simulations of the magnetic saturation remanence show that M has 
decreased only from 0.95 to 0.70. a m corresponds to an angular 
spread of about 45 degrees. Results of the simulations for intermedi- 
ate values of a g are summarized in Figs. 5d,e, revealing a direct 
proportionality between geometrical and magnetic disorder (a g , 
a m ) and a slow, close to parabolic, decay of M with geometric dis- 
order rj g . The persistent ferromagnetism in the quasi-2D structures 
appears contrary to the spin-glass behavior often expected in disor- 
dered nanoparticle- systems. 

Discussion 

For the arrangements of 15 nm Co particles studied here, our results 
show that short-range magnetic order with small domains dominates 
in the initial states, with the local magnetic order (ferromagnetic vs. 
antiferromagnetic) depending on the particle lattice (triangular vs. 
square). In contrast, at remanence after saturation, overall dipolar 
ferromagnetic order is extremely persistent even in case of a non- 
triangular lattice. We interpret our results as supporting the existence 
of amorphous dipolar ferromagnets: i.e., dipolar ferromagnetism 
exists in ID to quasi-2D structures even in the absence of underlying 



crystallinity. These observations are contrary to the magnetic glass 
behavior found in random 3D particle arrangements 1317 " 19 and are of 
direct relevance for the use of quasi-2D nanoparticle assemblies in 
magnetic devices such as thin elongated memories, tunnel junctions, 
and sensors. 

Methods 

Spherical s-Co nanoparticles, —15 nm in diameter and coated with oleic acid, were 
synthesized using a method adapted from Ref. 27, which involves thermal decom- 
position of a cobalt precursor in a mixture of hot organic solvent and surfactants at 
high temperature and under argon atmosphere using a Shlenkline setup to avoid 
oxidation. Cobalt carbonyl (Co 2 (CO) 8 , moistened with 5-10% hexane), trioctyl- 
phosphine oxide (TOPO, 99%), oleic acid (OA, 99%) and 1,2-dichlorobenzene 
anhydrous (DCB, 99%), were purchased from Sigma- Aldrich and used as received 
without further purification. All chemicals were stored and prepared for synthesis in a 
UNIlab MBraun glovebox. The cobalt precursor solution (0.54 g of Co 2 (CO) 8 
dissolved in 3 ml of DCB anhydrous) was injected rapidly into a refluxing DCB 
anhydrous solution (181°C) containing the surfactant mixture (15 ml of DCB in the 
presence of 0.2 g of OA and 0.15 g of TOPO). The solution was refluxed at 181°C for 
15 min under continuous stirring and then gradually cooled down before being 
transferred to an argon filled vial. One drop of the cobalt nanoparticle solution was 
then deposited onto a holey carbon-coated copper TEM grid and left to dry inside a 
glove box. 

Off- axis electron holograms were acquired at 300 kV in Lorentz mode in a Titan 
80-300 ST TEM with a charge-coupled device (CCD) camera and an electron biprism 
operated typically at 200 V. The experiments were performed at room temperature, 
and the acquisition time was 6-12 s. Initial magnetic states were recorded after 
loading the as-prepared samples into the TEM in almost field-free (below 6 mT) 
conditions. The direction of magnetization in each chain was subsequently reversed 
in-situ in the TEM by tilting the sample by ± 60° and turning on the conventional 
microscope objective lens to apply a magnetic field of 2 T to the sample, parallel to the 
direction of the electron beam. The objective lens was then turned off and the sample 
tilted back to 0° for hologram acquisition in field-free conditions and with the 
particles at remanence. Following this procedure, the different holograms were 
recorded with the chains magnetized in opposite directions. The mean inner potential 
was separated from the magnetic potential as described in [25]. 

The magnitude and orientation of each moment was measured by integrating the 
phase gradient components within each nanoparticle in each chain. For uniformly 
magnetized spherical particles, the integrated phase gradient lies in the same direction 
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and is directly proportional to the magnetic moment, with a factor of proportionality 
equal to n\i 0 /(2§ 0 ) = 0.954 rad'A^-nm" 1 , where |i 0 is the vacuum permeability, and 
4> 0 = 2.07* 10~ 15 T*m 2 is a flux quantum. The measurement is generally perturbed by 
the presence of neighboring particles proportionally to the magnitude and orientation 
of the field generated by surrounding particles. However, if the system is at equilib- 
rium, no torque is acting on any of the moments so that the local field and each 
moment are aligned; as a consequence, we do not expect significant perturbations in 
the measured orientation. 

The magnetic relaxation of the chains was simulated by an iterative approach based 
on a simplified version of the LLG equations, in which precession of moments around 
the local magnetic field was ignored (high damping coefficient). A simulation pro- 
ceeds in steps: 1) the simulated structure is created by assigning coordinates and 
diameters of the nanoparticles from an experimental micrograph; 2) each particle is 
assumed to be spherical, and with a magnetic moment proportional to its volume and 
with no magnetic anisotropy; 3) the initial orientation is chosen as either "saturated " 
where all moments are aligned (with a small random component), or "random" where 
all moments are pointing randomly; 4) the local magnetic field acting on the i-th 
particle is calculated as superposition of the magnetic fields generated by the 
remaining N-l particles; 5) the i-th moment is rotated towards the local field by an 
amount proportional to the torque, i.e. to the cross product between moment and 
local field; 6) a new moment topography is obtained, and the procedure is repeated 
until the total residual absolute torque becomes smaller than a pre-defined threshold. 
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